A combined theoretical and experimental study of the low temperature properties of 

BaZr0 3 
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Low temperature properties of BaZrC>3 are revealed by combining experimental techniques (X-ray 
diffraction, neutron scattering and dielectric measurements) with theoretical first-principles-based 
methods (total energy and linear response calculations within density functional theory, and effective 
Hamiltonian approaches incorporating/neglecting zero-point phonon vibrations) . Unlike most of the 
perovskite systems, BaZr03 does not undergo any (long-range-order) structural phase transition and 
thus remains cubic and paraelectric down to 2 K, even when neglecting zero-point phonon vibrations. 
On the other hand, these latter pure quantum effects lead to a negligible thermal dependency of the 
cubic lattice parameter below ~ 40 K. They also affect the dielectricity of BaZrC>3 by inducing an 
overall saturation of the real part of the dielectric response, for temperatures below ~ 40 K. Two fine 
structures in the real part, as well as in the imaginary part, of dielectric response are further observed 
around 50-65 K and 15 K, respectively. Microscopic origins (e.g., unavoidable defects and oxygen 
octahedra rotation occurring at a local scale) of such anomalies are suggested. Finally, possible 
reasons for the facts that some of these dielectric anomalies have not been previously reported in 
the better studied KTaC>3 and SrTiOs incipient ferroelectrics are also discussed. 
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I. INTRODUCTION 

BaZrOa is a ceramic oxide of the perovskite fam- 
ily structure with a large lattice constant, high melt- 
ing point, small thermal expansion coefficient, low di- 

luctivity (see, e.g. 
and references therein). 
These afore-mentioned properties make BaZr0 3 (i) a 
very good candidate to be used as an inert crucible 
in crystal growth techniques 7 ' 11 , (ii) an excellent ma- 
terial for wireless communication a 1 1 10 and (iii) a very 
good substrate in thin film deposition^. BaZrC>3 is 
also one of the two parent compounds of the (Pb-free 
and thus environmental-friendly) Ba(Zr,Ti)03 solid so- 
lutions, which is promising for manufacturing high Q 
materials with a variety of applications in microwave in- 
dustrjii 3 .. 

Interestingly, properties of BaZrOa have been mea- 
sured as long as 40 years ago, as well as very recently 
iSii, but only at room or higher temperature (to the 
best of our knowledge). Similarly, we are not aware 
of any calculation (either from phenomenological theory 
or first-principles calculations) predicting the dielectric 
properties of BaZrC>3. As a result, low-temperature di- 
electric properties of BaZrOs have never been investi- 
gated, despite the fact that many unusual effects are 
known to occur in some perovskite materials between 
and 50 K. One drastic example of such effects is the 
(temperature-independent) plateau and large values of 
the real part of the dielectric response in KTa03 and 
SrTiC>3, which arise from the quantum- induced suppres- 
sion of ferroelectricity in these materials. Other exam- 
ples are the anomalous peaks observed around 10-50 K for 
the imaginary part of the dielectric response in KTa03, 



K(Ta,Nb)0 3 , (Pb,La)Ti0 3 :Cui^ia, SrTiOjfiSA that are 
neither associated with structural phase transition nor do 
have a corresponding peak in the real part of the dielec- 
tric response (which conflicts with the well-established 
Kramers-Kronig relations!!). 

The aim of this article is to investigate the low- 
temperature behavior of BaZr0 3 from measurements and 
first-principles-based simulations. We report several un- 
usual features in the real and imaginary parts of the di- 
electric responses while no long-range ferroelectric, anti- 
ferroelectric or antiferrodistortive structural phase tran- 
sition occurs in BaZrOs down to 2K. Discussions and 
similarities/differences between BaZrOs and the (better- 
studied) KTaOs and SrTi03-related materials are also 
indicated to better understand the low-temperature di- 
electric anomalies reported in several perovskites. 

This article is organized as follows. Sec. ITT1 describes 
the experimental and theoretical methods we used to in- 
vestigate BaZrOs. Sec. IIIII reports the measurements 
and predictions for structural and dielectric properties. 
Finally, Sections IIVI and provide a discussion and 
conclusion, respectively. 



II. METHODOLOGY 

A. EXPERIMENTAL PROCEDURES 

Powder samples of BaZrOs were synthesized by solid 
state reaction by calcination at 1400K and sintering at 
1600K starting from stoichiometric amounts of the corre- 
sponding oxides (BaCOs, Zr02). The synthesized sam- 
ples were well crystallized and no presence of parasitic 
phases was evidenced by X-ray diffraction (XRD) and 



2 



TABLE I: The LDA-derived H c ff parameters in atomic units for BaZrC>3 following the notation in Ref [29]. 
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chemical analysis. The temperature dependence of the 
dielectric constant was measured at various frequencies in 
a temperature range from 5 K to 300 K using a Hewlett- 
Packard 4192A impedance analyzer and a cryostat with 
an estimated precision of 0.1 K. These measurements 
were performed on ceramic samples which were polished 
and cleaned, and sputtered gold electrodes were applied. 
Samples were annealed at 800 K and slowly cooled in or- 
der to eliminate strains caused by polishing. Powdered 
samples were used for the diffraction experiments. X-ray- 
diffraction measurements were performed on a high ac- 
curacy, two-axis diffractometer in a Bragg-Brentano ge- 
ometry using Cu-Kb wavelength issued from a 18-kW 
rotating anode generator, with diffraction angles preci- 
sion better than 0.002 deg. The neutron powder diffrac- 
tion patterns were collected at temperatures between 300 
and 2 K on the 3T2 high resolution goniometer on a ther- 
mal source (1.227 A) using the Orphe reactor facilities at 
Laboratoire Leon Brillouin (Saclay, France). Structural 
Rietveld refinements on both X-ray and neutron patterns 
were carried out with the XND software. 



B. THEORETICAL APPROACHES 

In this study, two different (direct) first-principles 
codes, as well as a first-principles-derived technique, were 
used to obtain and/or extract various information. One 
of the two first-principles codes is denoted ABINIT 18 . 
We took advantage of its implementation of the linear re- 
sponse theory to compute the phonon dispersion curves 
of BaZrC>3 at its low-temperature experimental lattice 
constant. Some technical details are as following. In 
this method we used Teter extended norm-conserving 
pseudopotentialsiSi and the local-density approximation 
(LDA) 2 ». The exchange-correlation functional was ap- 
proximated using Perdew-Zunger parametrization^i of 
Ceperley- Alder dataS. The Ba 5s 2 , Ba 5p 6 , Ba 6s 2 , Zr 
4s 2 , Zr 4p 6 , Zr 4d 2 , Zr 5s 2 , O 2s 2 , and O 2p 6 electrons 
are treated as valence electrons. A plane-wave cut off 
of 100 Ry was used to obtain convergence of the ground 
state total energy. Phonon frequencies were found to 
converge for dynamical matrices calculated on a 6x6x6 
Monkhorst-Pack grid 2 ^. 

The second first-principles code used in our study is de- 



noted CUSBSi. It also implements the LDASi, and the 
Ceperley- Alder exchange and correlation 22 as parameter- 
ized by Perdew and Zunger— . We also used the same va- 
lence electrons (indicated above) as in the ABINIT code. 
On the other hand, the pseudopotentials are those given 
by the Vanderbilt ultrasoft scheme 2 ^,, and the plane- wave 
cutoff is chosen to be 25 Ry, which leads to converged re- 
sults of physical properties of interest. The CUSP code 
is used to calculate the energetics of BaZr03 that are as- 
sociated with the rotation of the oxygen octahedra, at its 
low-temperature experimental lattice constant. It is also 
used to derive, at this specific lattice constant, the 18 
parameters of the first-principles-based effective Hamil- 
tonian (i? e ff) approach developed in Ref. |2(J These pa- 
rameters are given in Table I. 

We also perform simulations using such H c tf technique 
(and its parameters) to go beyond the abilities of direct 
first-principles techniques, namely to investigate finite- 
temperature properties of large BaZr03 supercells. The 
degrees of freedom of this Hamiltonian are the so-called 
local modes (that are directly proportional to the spon- 
taneous electrical polarization), and the strain variables 
(that characterize the crystallographic phase). The to- 
tal energy, E tot of H c ff contains five different interac- 
tions between local modes and/or strains, namely, a 
local-mode self energy, a long-range dipole-dipole inter- 
action, a short-range interaction between local modes, 
an elastic energy and an interaction between local-modes 
and strain — . Monte-Carlo (MC) simulations are per- 
formed using E to t in two different schemes: classical 
Monte Carlo (CMC)2£, which neglects zero-point phonon 
vibrations and, path-integral quantum Monte Carlo (PI- 
QMC)^2i22^£, which includes these quantum-mechanical 
zero-point motions. Consequently, comparison between 
the results of these two different Monte-Carlo techniques 
will allow the determination of quantum effects on struc- 
tural and dielectric properties of BaZrOa. 12x12x12 su- 
percells (corresponding to 8,640 atoms) are used in all 
Monte-Carlo simulations. We typically perform 30,000 
MC sweeps to thermalize the system and 70,000 more 
to compute averages, except at low temperatures in PI- 
QMC where more statistics is needed. 

In PI-QMC, each 5-atom cell interacts with its images 
at neighboring imaginary times through a spring-like po- 
tential (representing the zero-point phonon vibrations as 
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implemented in PI-QMC formulations), while all the 5- 
atom cells interact with each other at the same imaginary 
time through the internal potential associated with E to t ■ 
The product TP, where T is the simulated temperature 
and P is the number of imaginary time slices, controls 
the accuracy of the PI-QMC calculation^!. In all our 
simulations we use TP=600, which we find leads to suffi- 
ciently converged results. Outputs of the PI-QMC simu- 
lations thus contain local- modes (t) , where i runs over 
the 5-atom unit cells of the studied supercell while the 
imaginary time t ranges between 1 and P. (Note that 
t = P = 1 corresponds to CMC simulations.) Strain 
variables are another output of MC simulations. 



III. RESULTS 
A. STRUCTURAL PROPERTIES 

XRD and neutron scattering indicate that BaZr03 is 
cubic (and thus paraelectric) down to 2 Kelvin (the min- 
imal temperature accessed during our neutron scatter- 
ing). Note that we are not aware of any previous mea- 
surement investigating the low-temperature, rather than 
high-temperatur o 2 i 3 i 32 , properties of BaZr03. 

Our experimental result is consistent with the fact that 
we further find that our PI-QMC simulations using the 
effective Hamiltonian approach provide a vanishing su- 
percell average of the local-modes down to OK and a 
resulting cubic ground-state. Interestingly, our CMC 
simulations also predict such paraelectric ground-state. 
This indicates that, unlike KTaOg^ 3 - 3 ! and SrTiOa±£i22 
for which quantum effects suppress ferroelectricity in fa- 
vor of paraelectricity, zero-point phonon vibrations do 
not affect the symmetry of the ground-state in BaZr03. 
On the other hand, Fig.^ - which displays the predicted 
and measured temperature evolutions of the cubic lattice 
constant of BaZr03 - clearly shows, via the comparison 
of CMC and PI-QMC results, how zero-point phonon 
vibrations quantitatively affect structural properties of 
BaZr03. Below 100 K, the zero-point phonon vibrations 
tend to prevent the lattice constant from decreasing when 
decreasing the temperature. Quantum effects thus in- 
crease the lattice constant with respect to classical pre- 
dictions, with this increase becoming more pronounced as 
the temperature decreases. As a matter of fact, PI-QMC 
calculations result in (1) a dramatic change of thermal 
expansion below versus above 100 K, (2) a lattice con- 
stant that is nearly temperature-independent below ~ 
40 K, and (3) a lattice parameter that is larger by 10~ 3 
A from its CMC result at the lowest temperatures. Note 
that item (1) is clearly confirmed by our neutron scatter- 
ing data (reported in Fig. [2 for 300, 100 and 2K) that 
shows a relatively rapid ~ 10~ 3 A decrease of the lattice 
constant between 300 K and 100 K, while the difference 
in lattice parameters between 100 and 2 K is as small as 
3 x 10 _4 A. Furthermore, items (1) and (2) have also pre- 
viously been observed in the incipient KTa03 system^ 5 . 
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FIG. 1: Evolution of the cubic lattice constant as a function 
of temperature in BaZrOs. Solid squares and solid circles cor- 
respond to the CMC and PI-QMC predictions, respectively, 
while the open symbols display the results of neutron scat- 
tering. Note that the CMC predictions exactly agree with 
the experimental data at low temperature as a consequence 
of our choice of the ao parameter for the effective Hamiltonian 
approach (see Table I). 



(but for different critical temperatures). Moreover, item 
(3) provides a measure of the (relatively small) quantita- 
tive effects of zero-point phonon vibrations on the lattice 
parameter per se (rather than on the thermal expansion) 
of BaZr0 3 . 

Fig. d also indicates that above 100 K, one can safely 
use the approximation that the cubic lattice param- 
eter linearly depends on temperature with the ther- 
mal expansion coefficient being 0.27 xl0~ 5 (1/K), 0.23 
xlO~ 5 (l/K) and 0.47 xlO- 5 (l/K) for the CMC, PI- 
QMC and neutron scattering data, respectively. Our 
experimental value of 0.47 xlO~ 5 (l/K) compares rather 
well with the linear thermal expansion coefficient of 0.69 
x 10 _5 (1/K) previously measured by X-ray diffraction for 
temperature ranging between 273 and 873 K— . On the 
other hand, the thermal expansion of 0.27 xl0~ 5 (1/K) 
and 0.23 xlO _5 (l/K) predicted by our simulations un- 
derestimate the experimental values by a ratio of 2 — 3. 
This discrepancy between simulations and the measure- 
ments arises from the fact that the effective Hamilto- 
nian approach only incorporates the ferroelectric-related 
vibrations among the optical modes, while an accurate 
description of thermal expansion requires to take into 
account all phonon modes 3 ^. 

The effective Hamiltonian approach used in our study 
thus also neglects some phonon modes that may con- 
dense in BaZrOs, such as the R25 modes that are as- 
sociated with the rotation of the oxygen octahedra — . 
To check such possibility, we decided to compute the 
whole phonon dispersion of cubic BaZr0 3 , using the lin- 
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FIG. 2: Phonon dispersion of cubic BaZrOs, from first prin- 
ciples calculations. 
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FIG. 3: Oxygen octahedra rotational energetics (per 5 atoms) 
versus the angle of the rotation about the [001] direction from 
first principles calculations. The zero energy correspond to 
the paraelectric state. Solid line is a guide for the eyes. 



FIG. 4: Theoretical (using CMC and PI-QMC simulations) 
and experimental (for different frequencies) real part of the 
dielectric constant. Solid lines show the fit of the data into 
the Barret relation. The arrows indicate the position of the 
two observed dielectric anomalies that are discussed in the 
text 



pies CUSP program^) versus the angle of the octahedra 
rotation about the [001] direction: the minimum energy 
(occurring for an angle about 4 deg) is only ~ 1.5 meV 
deeper than the energy of the paraelectric phase (associ- 
ated with a zero angle). Such difference in energy corre- 
sponds to a rather small temperature of wl7K. It is thus 
highly possible that zero-point phonon vibrations pre- 
vent the occurrence (down to 2 K) of the macroscopic cu- 
bic paraelectric-to-antiferrodistortive phase transition, 
which would explain why our low-temperature XRD and 
neutron scattering experiments do not reveal any addi- 
tional peak related to a doubling of the unit cell. This 
quantum-induced suppression may act either globally 
(i.e., different unit cells do exhibit some rotation of their 
oxygen octahedra, but these rotations are not long-range 
correlated) or locally (i.e., there is no rotation of the oxy- 
gen octahedra in any unit cell). We will come back to this 
point, and to this «17K temperature, when discussing 
the results on dielectric properties. 



ear response theory as implemented in the first princi- 
ples ABINIT codei£. The results are shown in Fig. 
and indeed confirm the previous ab-initio prediction^ of 
the condensation of the R-point zone-boundary mode. 
(Note that these phonon calculations also confirm that 
BaZrOs does not exhibit any ferroelectric instability at 
the r point, even when quantum statistics are neglected). 
One has to realize, though, that the instability associated 
with these antiferrodistortive motions is rather weak, as 
demonstrated by the relatively small negative value of 
the frequency. This is further evidenced in Fig. |3 that 
displays the total energy (as predicted by the first princi- 



B. REAL PART OF THE DIELECTRIC 
RESPONSE 

Fig. ^ shows our experimental determination of the 
real part of the dielectric response, si, versus tempera- 
ture for different frequencies. The real part of the di- 
electric response of BaZrOs increases as the temperature 
decreases from 300 K for any given frequency. It then ex- 
hibits an overall low-temperature saturation to a plateau 
having a value that is slightly dependent on the frequency 
(ranging between 55 and 58 for frequency varying be- 
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TABLE II: Coefficients of the fit of the real part of the dielectric response into the Barrett relation. All the parameters of 
all the fits are allowed to relax, except the parameters C and T s for the fit of the PI-QMC data that are kept frozen to their 
deduced values from our experiments at the lowest frequency. Such freezing is necessary to allow the fit to converge because of 
the relatively large fluctuation inherent to PI-QMC simulations at low temperature. 

f[kHz] 5 10 40 100 200 500 800 1000 PI-QMC 

C 2900±185 2900 ±170 3200±200 3400±240 3600±280 4200±370 4600±440 4800±480 2900 

T S [K] 81±2 80±2 80±2 80±2 80±2 77±3 76±3 75±3 81 

T [K] -99±10 -102±9 -114±10 -124±12 -135±13 -160±17 -175±19 -182±20 -114±21 

B 38.8±0.3 39.6±0.2 40.0±0.3 39.8±0.4 39.5±0.4 39.0±0.5 39.0±0.6 39.0±0.6 36.7± 1.4 



tween 5kHZ and 1MHz). The existence of such plateau 
has already been reported in the incipient ferroelectrics 
KTaOg^M and SrTiQj 1 ^, and is usually thought to be 
associated with zero-point phonon vibrations. To check 
this fact, we also report in Fig. 0] the predictions from 
the H e g approach using both CMC and PI-QMC tech- 
niques. It is obvious that using CMC, the real part of 
the dielectric response continuously increases to higher 
value when the temperature decreases while on the other 
hand the dielectric constant computed within PI-QMC 
tends to saturate, in overall, at low-temperatureiS. Our 
simulations thus prove that quantum fluctuations play 
an important role on the low-temperature dielectric re- 
sponses of BaZr0 3 , as in KTaOj^i and SrTiOgiS^. 
However and as already mentioned in Section III. A, this 
quantum-induced modification of dielectricity in BaZr0 3 
is not accompanied by a suppression of ferroelectricity 
since BaZrOs does not have any ferroelectric instability 
even in the classical regime - as indicated by the fact 
that our CMC simulations lead to a cubic paraelectric 
ground-state and a lack of peak in e\. This distinguishes 
BaZr0 3 from both KTa0 3 and SrTi0 3 . Figure H also 
reveals that our PI-QMC simulations yield predictions 
that are slightly smaller in magnitude than our experi- 
mental results. In fact, this has to be expected since (i) 
the measurements displayed in Fig. 0] demonstrate that 
decreasing frequency leads to a decrease of the dielec- 
tric response at any temperature and (ii) the calculations 
correspond to the static regime, i.e. to a zero frequency. 
We can thus conclude that the PI-QMC results are in 
overall quite accurate, especially when realizing that di- 
electric coefficients are related to the derivative of the 
polarization (i.e., they are much more difficult to predict 
than properties that are directly proportional to polar- 
ization) and that the magnitude of e\ is rather small in 
BaZr0 3 (e.g., the low-temperature plateau has a value 
that is smaller by ~ 2-3 orders of magnitude than the 
corresponding ones in KTaQ 3 14 i 33 i 40 and SrTiOSiSi 41 .) 

To further analyze our results, we fitted them using 
the Barret relational 



C 



T s coth(-ff) - T 



B 



(1) 



is, the Curie temperature if quantum effects would not 
exist) and B is a constant independent of temperature. 
The resulting fits are indicated by means of solid lines in 
Fig. 21 and the coefficients of these fits are given in Ta- 
ble^ In addition to the fact that the quality of these fits 
is in overall rather good (thus, confirming, as in KTa0 3 
and SrTi0 3 , the relevance of the empirical Barret re- 
lation for describing quantum effects on dielectricity of 
perovskites), four features are particularly worth noticing 
related to these fits. First of all, To is (strongly) negative 
which confirms our theoretical findings that BaZr0 3 , un- 
like KTa0 3 and SrTi0 3 , does not exhibit any ferroelectric 
instability even in the classical regime. Secondly, the T s 
deduced from the experimental data is around 75-81 K, 
which is located in the temperature region for which PI- 
QMC predictions begin to significantly differ from CMC 
results for dielectric as well as structural properties (see 
Fig. 0| and Fig.^l. In other words, our simulations con- 
firm the physical meaning usually associated with T s , 
that is the temperature below which quantum effects 
play a non- negligible role on physical properties j&iM 3 .. 
Thirdly, and also unlike in KTa0 3 43 and SrTiO;A the 
B parameter can not be neglected to get good fits in 
BaZr0 3 . This latter difference is due to the fact that 
the low-temperature plateau is much larger in KTa0 3 
(around 4000) 14 i 33 i 4 9 and SrTi0 3 (around 20,000)iMI 
than in BaZr0 3 (around 55, see Fig. 0}, or equivalently 
that T is positive in KTaOg 3 ^ and SrTiOsii while be- 
ing negative in BaZr0 3 . Finally, two fine structures, ex- 
isting in the experimental data, deviates from (and seem 
to be superimposed with respect to) the Barrett fit. More 
precisely, one hump appears in E\ around 50 K at the low- 
est used frequency versus 65 K for the highest frequency, 
while a second hump shows up around 15 K for (more-or- 
less) any frequency. The magnitude of this second hump 
decreases when decreasing frequency. Note that these 
relatively small humps do not appear within our CMC 
simulations, and that the large fluctuations inherent to 
the PI-QMC approach^ do not allow us to assert if these 
humps are also predicted from our quantum simulations. 



where C and T s are the so-called Curie constant and 
saturation temperature, respectively Tq is in- 

terpreted as being the classical Curie temperature (that 



6 



0.012 



10 kHz - 
40 kHz 
100 kHz 
200 kHz 
500 kHz 
800 kHz 
1 MHz 




60 120 180 240 300 
T(K) 



FIG. 5: Measured (a) loss angle, tan<5, and (b) imaginary part 
of the dielectric response, £2 at different frequencies versus 
temperature. 



C. IMAGINARY PART OF THE DIELECTRIC 
RESPONSE AND DIELECTRIC LOSS 

Fig [S] displays the imaginary part of the dielectric re- 
sponse, £2 and the loss angle (as given by tan S = £2/^1) 
of BaZrOa versus temperature for frequencies ranging be- 
tween 10 kHz and 1MHz. All the results displayed in 
Fig are from measurements since simulating loss on 
an ab-initio level is one of the most interesting and diffi- 
cult challenges that remain to be accomplished nowadays. 
One can first notice from Fig [5] that both £2 and tan 5 
are very weak for the whole investigated frequency range, 
as also observed in Ba(Zn 1 / 3 Ta 2 /3)03-BaZr03 solid so- 
lutions 13 . BaZrC>3 can thus be a material of choice to 
design high-Q compounds. Furthermore, our measure- 
ments show the existence of two loss anomalies, namely a 
peak in £2 (and tan 6) around 50-65 K and a continuously 
increasing £2 (and tan<5) when decreasing temperature 
below 25 K. These two anomalies are correlated with the 
humps seen in £1 since they occur at similar temperatures 
and since they behave in a similar fashion. For instance, 
both the peaks in £1 and £2, occurring around 50 K at 
low frequency, shift to higher temperature when increas- 
ing the frequency. Similarly, the anomalies observed in s± 
and £2 at very low temperature get more pronounced at 
higher frequency. Interestingly, the well-known Kramers- 
Kronig relations 1 ? imply that any peak/anomaly in E\ 
should be accompanied with a peak/anomaly in £2 (and 
vice-versa) at the same temperature. Such correlation 



is thus indeed satisfied in BaZrC>3, but surprisingly, does 
not seem to hold for pure SrTi03 and KTa03. More pre- 
cisely, £2 has been found to exhibit a peak in these two 
latter materials for a temperature around 30-50 K, but 
no corresponding peak has been seen in the real part of 
the dielectric response. One possible reason for this lack 
of observation in pure KTa03 and pure SrTiOs may be 
due to the fact that these two materials (unlike BaZrOs) 
have a large overall real part of the dielectric respons o 16 i 33 
that "washes out" (i.e., prevents the observation of) weak 
superimposed peaks. 

Furthermore, the dielectric anomalies seen in BaZr03 
and reported in Figures 0] and El ar e also rather intrigu- 
ing since our XRD and neutron scattering do not detect 
any structural phase transition down to 2 K. A discussion 
about their possible causes is given in the next section. 



IV. DISCUSSIONS 

We believe that the enhancement in the real and imag- 
inary parts of the dielectric responses occurring at tem- 
perature around ~ 15 K is due to the activation of the 
oxygen octahedra rotation in BaZrOa. This belief is 
based on the fact that the minimum energy associated 
with antiferrodistortive motions is predicted to corre- 
spond to a temperature that is very close to the one at 
which these enhancements begin to occur, namely 17 K 
(see Fig. [3] and Sec. lIII All . The zero-point phonon vibra- 
tions may annihilate the long-range order of such rotation 
down to 2K, which would explain our X-ray diffraction 
and neutron scattering results. However, these quantum- 
induced effects likely can not prevent the octahedra ro- 
tation from occurring at a local scale (especially at very 
low temperature), which would be consistent with the 
fact that £2 increases when decreasing the temperature 
below 25 K (see Fig|5J|. 

Regarding the unusual dielectric features occur- 
ring around 50 K, it is important to realize that 
a frequency-dependent peak has previously been re- 
ported for tan S and £2 in several other perovskite sys- 
tems (e.g., BaTi0 3 :La, SrTi0 3 :La, SrTi0 3 , SrTi0 3 :Ca, 
K(Ta,Nb)0 3 , (Pb,La)Ti0 3 : Cu KTaOs) near a similar 
temperature (see, e.g. Refs Il4ll5llr>l34l and references 
therein). For instance, £2 peaks around 40 K for low fre- 
quency in KTaO^i^i. It is commonly believed that the 
reasons behind these peaks is the existence of unavoid- 
able impurity ions having a, different valence than the 
host atoms (see Refs. I44M5I and reference therein). In 
such a case, these peaks should become more pronounced 
when intentionally doping the sample with impurities, as 
consistent with the fact that adding up to 3% of Mn 2+ 
in KTa03 leads to a noticeable peak in the real part of 
the dielectric response around 40 K^. Observing and un- 
derstanding the effects of doping on physical properties 
of BaZr03 is thus of importance to confirm this (gen- 
eral) possibility, but goes beyond the scope of the present 
article. One particular previous studyi^ further stipu- 
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lates that it is a polaronic relaxation - that is, a coupling 
between the free charge carriers arising from the impu- 
rity ions and the lattice properties of the host material 
- that causes such weak dielectric anomalies. As done 
in Ref. ITU such possibility can be checked by extracting 
the temperature-dependency and frequency-dependency 
of relaxation time via an analysis of loss dynamics using, 
e.g, Cole-Cole formula lossii. Such analysis, to be accu- 
rate, requires the investigation of e\ and £2 under a range 
of frequency that is much wider than the one available 
for the present study. 

V. CONCLUSIONS 

In summary, we combined measurements with first- 
principles-based techniques to investigate the low- 
temperature properties of BaZr03. This system is found 
to be cubic and macroscopically paraelectric down to 2 K. 
Unlike the "better-studied" KTaOgS^a and SrTiOgi^ 
incipient ferroelectrics, the zero-point phonon vibrations 
do not suppress ferroclectricity in BaZrC>3. In other 
words, this latter material is also paraelectric in the clas- 
sical regime. On the other hand, quantum effects lead 
to the saturation of the cubic lattice parameter below ~ 
40 K. 

Despite having no long-range-order structural phase 
transition, BaZr03 exhibits the following striking dielec- 
tric features: (1) the £1 real part of the dielectric response 
saturates in overall at low temperature (namely, below ~ 
40 K) to a value ~ 55. Our PI-QMC simulations are in 
rather good agreement with experimental data, and show 
that such saturation is caused by zero-point phonon vi- 
brations; (2) the temperature behavior of £1 can be well 
fitted by the empirical Barrett relation when allowing the 
B parameter to differ from zero in this fit; (3) two peaks 
or fine structures are observed in £1 around 50-65 K and 
15 K, respectively. The first peak shifts to lower tem- 
peratures when decreasing the frequency, while the sec- 
ond one occurs at around the same temperature but de- 
creases in magnitude when decreasing the frequency; (4) 



these two fine structures are associated with anomalies in 
the imaginary part of the dielectric response, £2, which 
peaks at around 50-65 K while suddenly and continuously 
increases when decreasing temperature below 15 K. 

By comparing with previously reported data in other 
perovskites and adopting some related interpretations, 
we propose that the highest-temperature dielectric 
anomalies are caused by defects like oxygen vacancies 
and/or unavoidable impurity ions such as Fe 3+ — which, 
e.g., are the source of free charge carriers, that can inter- 
act with a soft lattice to create a polaronic state respon- 
sible for the dielectric anomaly^. Furthermore, our (di- 
rect) first-principles calculations suggest that the lowest- 
temperature dielectric anomalies may result from local 
rotation of the oxygen octahedra. 

In order to acquire a deeper knowledge of perovskites, 
it is interesting to compare items (1-4) with correspond- 
ing features in KTa03 and SrTi03. For instance, £1 also 
saturates at low temperature in KTa03 and SrTi03, but 
with a much higher value of the plateau (as a result of the 
quantum-induced suppression of ferroelectricity in these 
latter compounds). Such high value of the plateau ex- 
plains why the B coefficient is generally omitted in the 
Barrett fit of KTa03 and SrTi03. It is also highly plau- 
sible, as we believe it, that such high e\ prevents the 
observation of the hump that should be associated (ac- 
cording to the Kramers-Kronig relationii) with the peak 
of £ 2 seen around 30-50 K in KTaOgi^ and SrTiO;j±£J£. 

We hope that our work stimulates further investiga- 
tions aimed at checking our suggestions, in particular, 
and understanding dielectric anomalies in perovskites, in 
general. 
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